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The nuclear modification factor Raa of non-photonic electrons in Au + Au collisions at ^/s NN = 
200 GeV is studied by considering the decays of heavy-flavor hadrons produced in a quark coalescence 
model. Although an enhanced A c /D° ratio is predicted by the coalescence model, it is peaked at 
small transverse momenta (~ 2 GeV) due to the large difference between heavy and light quark 
masses. As a result, the enhanced A c /D° ratio, which is expected to suppress the electron Raa as 
the branching ratio of A c decay into electrons is smaller than that of D° , does not lead to additional 
suppression of the electron Raa at large transverse momenta (> 5 GeV), where the suppression is 
mainly due to heavy quark energy loss in produced quark-gluon plasma. Also, the enhanced A^/B 
ratio predicted by the coalescence model has even smaller effect on the non-photonic electron Raa 
as bottom baryons and mesons have similar branching ratios for semi-leptonic decays into electrons. 

PACS numbers: 25.75.-q, 25.75.Cj, 25.75.Dw 



Heavy-flavor hadrons are useful for probing the prop- 
erties of the dense matter produced in relativistic heavy- 
ion collisions. Although D mesons were measured in 
d + Au collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) [Ij, there is not yet direct measurement of 
open charm or open bottom hadrons in heavy-ion colli- 
sions. Instead, heavy-flavor hadron production in heavy 
ion collisions was studied through measurements of non- 
photonic electrons from their decays @. The transverse 
momentum (p T ) spectrum of these electrons in Au + Au 
collisions at ^/s NN = 200 GeV were measured for sev- 
eral centralities by the PHENIX Collaboration . The 
resulting nuclear modification factor Raa, defined by 
the ratio of the p T spectrum of non-photonic electrons 
in heavy ion collisions dN AA /dp T to that from proton- 
proton collisions dNp p /dp T multiplied by the initial num- 
ber of NN binary collisions (N^), i.e., 



Raa 



dN AA /dp 7 



(1) 



shows that the production of heavy-flavor hadrons in 
heavy-ion collisions is suppressed as much as that of pi- 
ons. This has raised a challenging question on heavy 
quark interactions in medium as perturbative Quantum 
Chromodynamics (pQCD) predicts that heavy quarks 
lose less energy in the quark-gluon plasma than glu- 
ons and light quarks, which is responsible for the ob- 
served large suppression of high p T pions. Various ideas 
have been suggested to explain the observed small non- 
photonic electron Raa, and these include the introduc- 
tion of very large medium opacity Q , two-body [1, d, 0| 
and three-body collisional energy loss y|, and the colli- 
sional dissociation of heavy mesons in medium Q . It is 
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also proposed that an enhancement of the A c /D° ratio 
would reduce the Raa of non-photonic electrons [l(| EH • 
This is based on the observation that the branching ratio 
for A c — > e decay is smaller than that for D° — > e de- 
cay. Therefore, if the A c /D° ratio is enhanced as in the 
case of p/n and A/K° ratios, the Raa of non-photonic 
electrons would be smaller than that only due to heavy 
quark energy loss. In Ref. [ll|, it was claimed that this 
effect may explain about 20% of the suppression of Raa- 
Since the origin of the A c / D° enhancement is unknown, 
this ratio is assumed in Refs. [l(| HH to have either the 
same p T dependence as in the A/K data or a Gaussian 
form. 

Recently, the enhancement of A c yield in relativistic 
heavy ion collisions was suggested as a possible test of the 
diquark model for the structure of heavy baryons [l2T.[l3l|. 
Estimates based on the quark coalescence model for 
heavy-flavor hadron production show that the existence 
of diquarks in both the quark-gluon plasma and the A c 
would enhance the A c /D° ratio by a factor of ten and 
five, respectively, compared to those of PYTHIA simu- 
lation for pp collisions and of a simple thermal model. 
(See, e.g., Ref. [lij for the quark coalescence model.) 
Furthermore, it was shown that the enhancement of the 
A c /D° ratio could be better seen at low p T region where 
heavy quark fragmentation is less important. Therefore, 
the enhancement of A c /D° ratio is expected to result in 
the suppression of Raa only at low p T re gion . This is 
in contrast to the assumption of Refs. fiol llll| that the 
A c /D° enhancement continues to the large p T region. In 
this report, we calculate the non-photonic electron Raa 
based on the results from Ref. [la ]. We also discuss the 
contribution from the enhanced A^/B ratio to the non- 
photonic electron Raa- 

For calculating the electron spectrum from heavy 
hadron decays, information on the branching ratios of 
scmi-leptonic decays of heavy hadrons into electrons is 
needed. The experimental data complied by the Particle 
Data Group [15| for the branching ratios of these de- 
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Decay channel 




D° 


D+ 


D+ 


A c 


BR(e + + anything) 


6.53 


±0.17% 


16.0 ±0.4% 


_|_g 

8 -5 % 


4.5 ±1.7% 


BR{Ke+u e ) 


3.58 


± 0.06 % 


8.6 ±0.5% 






BR(K*e + v e ) 


2.18 


±0.16% 


3.66 ± 0.21 % 






BR(r?£ + ^) 








2.9 ±0.6% 




ER(r)'t+u £ ) 








1.02 ±0.33% 




BR(Ae+i/ e ) 










2.1 ±0.6% 


BR(pe + + anything) 










1.8 ±0.9% 



TABLE I: Branching ratios of major semi-leptonic charm hadron decays into electrons as listed by the Particle Data Group |l5| |. 



cays are, however, ambiguous as shown in Table [U One 
can estimate the effect of the A c /D° enhancement on the 
electron R AA roughly by considering semi-leptonic decay 
of heavy hadrons in each channel as in Ref . [ll[ , 



tion given in Ref. 



Raa 



{n e (D°) + n e (D+) + n e (D s ) + n e (A c )} 



AA 



{n e (D°) ± n e (D+) ± n e {D s ) + n e (A c )}_ 



where n e {H) is the number of electron produced by the 
decay of hadron H, and the number of binary collisions 
is included in the numerator. Since the branching ratio 
of A c decay into electron is smaller than that of charmed 
mesons, the enhanced A c /D° ratio indeed suppresses the 
electron Raa if the momentum difference between elec- 
trons and the decayed charmed hadrons is neglected. In 
Ref. [llj], it was claimed that if the enhancement factor, 
C = (A c /D°) AA /(A C /D°) pp , is 12, the electron R AA is 
reduced by about 20%. More quantitative understand- 
ing of the effect of enhanced A c /D° ratio on the R AA of 
non-photonic electrons requires, however, a better under- 
standing of the branching ratios of heavy hadron decays, 
including the decays to four-body final states, as well 
as the production mechanism of heavy-flavor hadrons in 
both pp and AA collisions. The latter is addressed in the 
present study. 

We use the coalescence model calculation presented in 
Ref. [HI], which is based on the coalescence formulas 



dN E 
dp B 



9 b- 



(2^) 6 (a 1 a 2 



V 2 



dp 1 dp 2 dp : 



dNx dN 2 dN 3 
dpi dp 2 dp 3 



x exp (-fefof - fc^crf) 8(p B -Pi-p 2 - P 3 ), 

(3) 

for baryon production from three-quark coalescence and 



dN M 

d P M 



9m- 



V 



dp x dp, 



exp (- 



-k 2 a 2 dq ) 5(p B 



dNx dN 2 
dp l dp 2 

Pi -P 2 ), 



(4) 



for both baryon and meson production from diquark- 
quark coalescence and from quark-anti-quark coales- 
cence, respectively. For the details on the model and 
parameters used in the calculations, we refer to Ref. [l3l |. 
For the heavy quark energy loss, we use the parametriza- 



L c = 0.8exp(-p/1.2)±0.6exp(-p/15), 
L b = 0.36 ± 0.84 exp(-p/10), 



(5) 



for charm and bottom quarks, respectively, where p is 
the transverse momentum in units of GeV. As shown in 
Ref. [HI , compared with charmed hadron production via 
fragmentation of charm quarks, the yield of D° mesons in 
heavy ion collisions is suppressed in the quark coalescence 
model as a result of enhanced production of A c . (See the 
small window in Fig. Q]for the A c /D° ratios.) This hap- 
pens in models based on both the three-quark coalescence 
and the diquark-quark coalescence. In Fig.[TJ we plot the 
electron spectrum from charm hadron decays for central 
Au + Au collisions at mid-rapidity (\y\ < 0.5). As a ref- 
erence, results from the decays of charmed hadrons pro- 
duced via 8- function fragmentation of charm quarks, i.e., 
charm quark and produced hadrons have same p T , with- 
out and with energy loss are given by the solid line and 
dashed line, respectively. These electron spectra include 
contributions from decays of D°, D + , D s , A c and their 
anti-particles, as production of electrons from charmed 
hadrons other than D° is not negligible, particularly if 
there is a large enhancement of the A c /D° ratio [lfj |. 
Compared with the case that these charmed hadrons are 
produced from fragmentation of charm quarks via a 8- 
function fragmentation, coalescence models are seen to 
give a somewhat suppressed electron spectrum. We note 
that the uncertainties in the branching ratios as shown 
in Table |I] can affect the electron spectrum. More discus- 
sions on these uncertainties can be found in Ref. [iH ]. 

Our results for the non-photonic electron R AA are pre- 
sented in Fig. [2] Here, we consider five cases for inves- 
tigating the role of the enhancement of the A c /D° ratio 
in the electron R AA . First, the solid line is the electron 
R AA for the case with charm quark energy loss in the 
(5-function fragmentation approximation as discussed in 
Ref. Q. Shown by the dotted line is the result from the 
coalescence model when only D mesons are allowed to 
be produced. If other charmed hadrons are allowed to be 
produced, the contribution from the D° mesons to the 
electron spectrum is suppressed and leads to the sup- 
pressed electron R AA as can be seen by the dot-dashed 
line. The dashed and dot-dot-dashed lines are obtained 



3 




FIG. 1: Spectra of electrons from decays of all charmed 
hadrons. Shown in the inset are the A c /D° ratios in frag- 
mentation and in both the three-quark and diquark-quark 
coalescence models. 



by allowing the decays of all charmed hadrons into elec- 
trons in diquark and three-quark coalescence, respec- 
tively, and are the main results of this study. Therefore, 
by comparing the dotted line and the dashed line, one 
can see that the enhanced A c / D° leads to suppression of 
the electron Raa- 

However, when we compare our result for the electron 
Raa with that of fragmentation with energy loss, the 
enhancement of the A c / D° ratio does not cause the sup- 
pression of the electron Raa- This is because the en- 
hancement of A c /D° ratio in the coalescence model ap- 
pears mainly in the region of p T < 5 GeV. In addition, we 
found that the crossover of the coalescence and fragmen- 
tation is at p T ~ 0.75 GeV. As a result, the suppression 
of the Raa for non-photonic electrons appears also at low 
p T region (p T < 3 GeV). Thus, the calculated Raa un- 
derestimates the data [3j at low p T and overestimates it 
at large p T . Including the decay of A c , whose branching 
ratio of decaying into an electron is not small in com- 
parison to that of D°, increases the electron Raa- We 
also found that the difference between the results of the 
three-quark and the diquark-quark coalescence models is 
not large. Therefore, in the coalescence model considered 
here, more quenched charm quark spectrum is needed to 
explain the small electron Raa observed in experiments. 

Our results that the enhancement of A c /D° ratio does 
not cause the suppression of Raa is in sharp contrast 
with those of Refs. [13, EH- The reason for this is due to 
the difference in the A c /D° ratio at large p T region. In 
both Refs. [T3, EH , a constant enhancement factor is as- 
sumed for the A c /D° ratio at large p T region. Besides the 
assumption that the A c /D° ratio is the same as the A/AT§ 
ratio, it is further assumed in Ref. that A C /D Q = 0.33 
for large p T , where no data for A/Kg exist. In Ref. plj . 
an even larger constant enhancement factor C ~ 12 was 
introduced to the A c /D° ratio that was assumed to have a 




P T (GeV) 

FIG. 2: The electron Raa in the central collision of Au+Au 
at ^/s NN = 200 GeV from charmed hadrons. The solid 
line is the result obtained from only fragmentation of heavy 
quarks with energy loss. The dot-dot-dashed (dashed) line is 
for the three-quark (diquark-quark) coalescence model, while 
the dot-dashed line is obtained with D° meson only in the 
diquark-quark coalescence model. The dotted line is from the 
toy model described in the text and the data are from Ref. @|. 



Gaussian shape in p T . This large enhancement of A c /D° 
at large p T would then lead to the suppression of Raa at 
intermediate and large p T (> 5 GeV). 

In our case, although the A c /D° ratio is enhanced com- 
pared to that in pp collisions, the enhancement is mainly 
at small p T and the A c /D° ratio remains comparable to 
that of pp collisions at large p T (> 7 ~ 8 GeV) region. 
Therefore, the effect of the A c /D° enhancement on the 
electron Raa cannot be seen in the intermediate electron 
p T region. Furthermore, charmed hadrons produced by 
coalescence with light quarks have larger values of p T 
than that of the charm quark inside a hadron, which re- 
sults in the shift of charmed hadron p T spectra to the 
larger p T region. Therefore, at large p T region, we have 
more charm hadrons than those produced by 5 function 
fragmentation of charm quarks. Although this effect is 
not large compared to the case for light hadrons, it in- 
creases the electron Raa- As a result, we have a larger 
Raa at intermediate electron p T compared to that from 
the fragmentation of quenched heavy quarks. 

As the electron p T increases, electrons from bottom 
hadron decays become important. The p T at which 
electrons from charmed hadron decays and from bottom 
hadron decays intercepts depends on the model for heavy 
quark energy loss in the quark-gluon plasma. Since h eavy 
quark energy loss is suppressed as its mass increases pj]] , 
this momentum becomes smaller once heavy quark en- 
ergy loss is included. Therefore, semi-leptonic decays of 
bottom hadrons should be considered at large electron 
p T region, and this would further increase the electron 
Raa- 

The decays of bottom hadrons into electrons are not 
well-known. There are two processes through which elec- 



4 



1.5 



0.5 









I I I I I I I 

fragmentation (energy loss) 




3-quark coalescence 




diquark-quark coalescence 


-— 


coalescence ( toy model) 








P T (GeV) 



FIG. 3: The electron Raa in the central collision of Au+Au 
at y / S]vjv — 200 GeV from bottom hadrons. Notations are 
the same as in Fig. [5] 

trons are produced from bottom hadron decays. One 
is the direct production of electrons from bottom hadron 
decays, such as B — > ev e X c , where X c denotes any charm 
meson. The produced X c can again decay into elec- 
trons. Therefore, bottom hadrons produce primary and 
secondary electrons. In addition, for vector mesons, B* 
and B*, they first decay into B7 and -B s 7, respectively, 
and electrons are then produced from the decay of re- 
sulting pseudoscalar mesons. Currently available exper- 
imental data show that the branching ratios of bottom 
hadrons decaying into electrons are nearly independent 
of the light quark content of bottom hadrons [15| , 

BR(B ± — > e + v e X c ) = 10.8 ±0.4%, 

BR(5° -> e + v e X c ) = 10.1 ±0.4%, 

BR(£ S — ► £ + v e D s + anything) = 7.9 ±2.4%, 

BR(A fc -> £-p e A+ + anything) = 9.9 ±2.6%. (6) 

The paucity of experimental data for the branching ra- 
tio of each channel in semi-leptonic decays of bottom 
hadrons makes it difficult to compute the electron spec- 
trum from bottom hadron decays. Nevertheless, one can 
estimate the effect of A& /B° enhancement to the electron 
Raa in a qualitative way by using the formula given in 
Eq. 12]). In Ref. [13j |. we have observed a large enhance- 
ment of the Ab /B° ratio in the coalescence model. How- 
ever, since the branching ratio of B — > e and A^ — > e 



are similar, enhancement of A;, baryon produces more 
electrons, which compensates almost the reduced num- 
ber of electrons from B decays. This can also be verified 
through Eq. (J5J). Namely, if we use the branching ratios 
given in Eq. ([6]), Raa changes from 1.0 with C = to 
0.95 with C — ► 00. This shows that the effect of the 
A(,/B enhancement on the electron Raa is very small 
compared to that for charm hadrons. This is shown ex- 
plicitly in Fig. [3l By considering the large uncertainty 
in the branching ratio of Aj, decays, the effects of A v /B° 
enhancement is small, in particular at large p T region, 
where the electron Raa is expected to be governed by 
bottom hadron decays. 

Although the multiplicities of multi-heavy-quark 
baryons are very small, production of heavy baryons con- 
taining one heavy quark and one or two strange quarks 
is not negligible. In the coalescence model of Ref. [13j |. 
we found that the sum of the multiplicities of Sq, Sq, 
and Hq is as large as 30 ~ 50% of the multiplicity of 
Aq, where Q stands for c or b. However, the branching 
ratios of these hadrons have not been measured and the 
electron spectrum from the decays of these baryons thus 
can not be estimated. 

In summary, we have studied the transverse momen- 
tum spectrum and the nuclear modification factor Raa 
of non-photonic electrons from heavy hadron decays in 
relativistic heavy ion collisions. Contrary to the models 
of Refs. (hIIIIIj which assume an enhanced A c /D° ratio 
at large p T region, the large enhancement of A c /D° ra- 
tio predicted in the coalescence model [l3} occurs mainly 
in the low p T region and the heavy quark fragmentation 
remains dominant at large p T region. As a result, no ad- 
ditional suppression to the electron Raa due to enhanced 
A c /D° ratio is obtained at large p T region. We have also 
estimated the role of Af,/B° enhancement in the elec- 
tron Raa and found that the enhancement of A^ baryon 
production does not affect the electron Raa produced 
by bottom hadron decays because of similar branching 
ratios of bottom meson and A& decays into electrons. 
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